Two complementary descriptions of intermittency 
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We describe two complementary formalisms designed for the description of probability density 
function (PDF) of the gradients of turbulent fields. The first approach, we call it adiabatic, describes 
PDF at the values much less than dispersion. The second, instanton, approach gives the tails of 
PDF at the values of the gradient much larger than dispersion. Together, both approaches give 
satisfactory description of gradient PDFs, as illustrated here by an example of a passive scalar 
advected by a one- dimensional compressible random flow. 
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Probably the most striking feature of developed tur- 
bulence is its intermittent spatial and temporal behav- 
ior. The structures that arise in a random flow manifest 
themselves as high peaks at random places and at ran- 
dom times. The intervals between them are character- 
ized by a low intensity and a large size. Rare high peaks 
are responsible for PDF tails while the regions of low in- 
tensity contribute PDF near zero. That physical picture 
prompts an attempt to describe intermittency at the level 
of a single-point PDF by two complementary approaches. 
The first approach was recently introduced to describe 
rare strong fluctuations as optimal fluctuations realizing 
probability extrema . Called an instanton approach, 
this formalism is based upon a path-integral representa- 
tion of conditional probability with optimal fluctuations 
being saddle points in the integral. A counterpart to the 
instanton approach is suggested here for the description 
of gradient PDF at small values, the approach is just an 
adiabatics when high-order spatial derivatives are consis- 
tently neglected. 

The center anomaly and tails of the gradient PDF are 
two sides of the same coin called intermittency which is 
the main target in modern turbulence studies. It this Let- 
ter, we demonstrate how both methods applied together 
give a consistent description of the gradient PDF. Note 
that the intermediate part of the PDF (which is beyond 
our approaches) where the matching of the asymptotics 
occurs is not that interesting because it is nonuniver- 
sal i.e. depends on the particular form of the pumping 
correlation function. The central peak and the tail are 
robust, their form provides the main information on the 
probability of both main body of the events and strong 
fluctuations. 

Let us show how such a description can be developed 
by using, probably, the simplest (yet nontrivial) turbu- 
lent problem of a passive scalar 9{x,t) advected by one- 
dimensional random flow u(x, t) which is smooth in space 
and white in time (this is a compressible version of a 
well-known Kraichnan model 0): 

dte + vve = KAe + (/>. (i) 



Both the velocity v and the source function are sup- 
posed homogeneous, Gaussian and (5-correlated in time: 
{(f){x,t) (l){x' ,t')) = x{\x — x'\)6{t — t') where x is some 
function that decays on a scale L, the value x(0) = P 
is the flux of 9^ . The correlation function of the veloc- 
ity may be defined by two parameters, typical velocity V 
and correlation length 

{v{x,t)v{x' ,t')) = [VL^-VL-^{x-x'f]5{t^t'). (2) 

At studying a simultaneous statistics, the coordinate- 
independent part drops out. We assume 3> L. 

Let us first implement a simple adiabatic approach ne- 
glecting diffusion term. Then for the single-point PDF 
V{uj,t) — {S[9x{x,t) — uj]) one obtains a closed Fokker- 
Planck equation 

r)V 8^V r)V 

where we denote T = x"(0) and D = VL^^ the vari- 
ances of 4>x and Vx respectively. That equation has an 
equilibrium steady solution 

V{uj) (X {T + Duj^)-^ (4) 

expected to be applicable for oj'^ <C P/k. Since T ~ P/L^ 
and the Peclet number Pe^ = DL^/4k is assumed to 
be large then (||) has a wide interval of validity. Note 
that T/D is a square gradient produced by the pumping 
during the typical stretching time D~^. For Pe » 1, 
T/D <^ P/k. Limiting solutions obtained at w"^ ^ T/D 
and at cj^ <C T/Z? by a time-separation procedure j|] 
coincide with (Q). 

Let us now describe the tail of the probability den- 
sity function V{uj) at uj^ ^ P/k. It is clear from (|^) 
that the correlation functions of the strain field a = 
are x-independent that is <j can be treated as a random 
function of time t only. To exploit that, it is convenient 
to pass into the comoving reference frame that is to the 
frame moving with the velocity of a Lagrangian particle 
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of the fluid 0,^. The Martin-Siggia-Rose action / for 
the n-th order moment of the gradient 6^ is H: 



1= / dtdxpdte- I Edt~ — \og{e^f (5) 



^ J dxidx2 Pix{xi2)p2 ' ^ (6) 

Assuming n ^ 1 we shall calculate the moment in the 
saddle-point approximation (0") — eyip^Iextr), see 
for the details. Here, lextr is to be calculated on the flow 
configuration (optimal fluctuation or instanton) that has 
to satisfy extremum equations for the action: 

dtO + axd^e - ndle = -i j dx'p{t, x')xix - x') (7) 
dtp + adxixp) + kOIp ^ -y5{t)d'{x) (8) 
a-D = 2iD / dxpxdxO (9) 



Here y — —in/6x{0,0)- For calculations, it is more 
convenient to use this auxiliary parameter instead of 
6x{0, 0). The boundary conditions are 6{x, — oo) = and 
p{x,+0) = 00]. The solution of (@,||) can be sought in 
the following form 

e^f{T{t),xy/^),P^^^g{T{t),xy/^^ (10) 
dtw = ~2aw, dtT = -w, w{0) = 1, t(0) = 0. (11) 



The functions / and g satisfy the following equations 



drg-Kdlg^y6{T)S'iO 



(12) 



drf + ndlf^^ rd^'g{r,axi^-^r4^] (13) 
The solution can be found in the Fourier-representation 



g{T, k) — ikye 
fir, k) = -ky 



— Kk^T 



(14) 



dr' 



=x(fcv/^)e'''='(--2-'' (15) 



o r dr' [°° dk , 
J TO Vw(t') J-oo27r 

XX (k^/My)) e^'"'^'"' (16) 

Here tq is the maximal value for r, which is determined 
by the moment when the following integral diverges 



t(r) = - 



^ dr 
w{t)' 



(17) 



In the following we will work in the dimensionless units. 
We put D = 1, P = x(0) = 1, and L = I. Then 



K = l/(4Pe^), where Pe is the Peclet number. We be- 
lieve Pe > 1. Calculating dJiO^Qi) from (|l|) we obtain 
the following self-consistency condition for y 



<t>{x) 



dr 



w^(t) \Pe'^w{T) 
dk 



27r 



k x(fc)exp 



k'^x 



(18) 
(19) 



The function 4'{x) has the following asymptotics: 4>{x) — > 
1 as a; ^ and 4>{x) ~ as a; — > oo if x(A = 0) 0- 

Note, that x{^) ^ 0, hence ip{x) is a monotonic decreas- 
ing function. One may keep in mind some particular form 
of x(a;), say exp(-cc2LV2). Then, (l){x) = {l + x)'^/^. 

Now we derive a closed equation which describes the 
evolution of w which is the square root of the solution 
inverse width. One can do that directly from ( [T^ ) sub- 
stituting there a = w'{t)/2 which is the consequence of 
(|TT|). One obtains an integral equation which is equiva- 
lent to some third-order ordinary differential equations. 
The order can be then reduced by one due to the conser- 
vation law (^). However, it is more instructive to derive 
the same equation on w in a different manner: since we 
are looking for the extremum of the action (|^) we can 
substitute there all the fields as the functionals of w and 
then make a variation with respect to w. We have 
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1 r° dr 




4^0 w 


(r) 





dr 



w' - 1^ 



(20) 



Varying with respect to w, we obtain 



2w 



2 

w 



Pe^vj^ 



Pe^w 



Pe^w 



= 



(21) 



Equation ( pl| ) can be rewritten as a Hamiltonian system 
with the momentum P — w' jw and the Hamiltonian 



P^w 4|y|' 



Pe\ 



2 

w 



(22) 



Initial conditions for (21) are w{Q) = 1 and ^'(O) 



-2(2n - 1). The latter is readily derived from ( p4n8|) . 
We should satisfy also a final condition. Indeed, (17 20| ) 
require 2 as t — > tq, otherwise the integral contri- 

bution to the action is infinite. Our aim is to find such 
a value of y that the solution of (|2l|) satisfies the rela- 
tion (p^). In other words, we can divide the task into 
two parts. First, to find solution of equation ( pT| ) with 
arbitrary y and the given boundary conditions. Second, 
given w, to solve the algebraic equation ( p^ ) which deter- 
mines y. Note that finite positive value of w' at r — > tq 
implies tq = oo. 
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Let us start the first part of our program, which is 
solving (21). Since ( p2|) explicitly depends on r, then H 
is not conserved and (|2l|) can not be solved explicitly for 
an arbitrary (j). In the limit n <C Pe^, it is possible nev- 
ertheless to describe the solution with enough details to 
recover it's dependence on the parameters n and Pe. 

The evolution of w can be divided into three parts. 
During the initial stage, when r is close to zero, w is of 
order unity. Therefore, Pe^w <C t and we can substitute 
4> in (^) by its asymptotic value 1. Thus, during that 
stage 7? is a constant which we denote Hq. One finds 
Ho = H{0) = 2(2n - l)^ - 2 + A\y\^. Since n > 1 and, 
as we shall see below y <C 1, we have Hq w 8n^. The 
equation for w can be then readily derived from (122 



'/ = -v/4 + 2i/ow-8|?;|2/w 



(23) 



Now let us consider the final stage. Since w' ^ 2 as 
T oo, we can write w « 2r. It gives T/{Pe^w) ~ 
l/(2Pe^) <C 1. Thus, like for the first stage, we can re- 
place by 1, and therefore the energy is a constant. It 
follows from ( p2| ) that during that stage Hw <^C 1. Hence, 
w satisfies the following equation 



w' ^ ^4-8|?;|2/w 



(24) 



Since w' is negative during the first stage and positive 
during the third one, then it has to turn into zero at 
some refiection time r*. Around that time, there should 
exists an intermediate part of evolution, which matches 
the two above asymptotics. We will see below, that this 
stage gives the main contribution into dl^). During that 
stage one has r/iPe'^w) > 1, so that H is not conserved 
but decreases from 2n^ to 0. It will be important for us 
that H is a decreasing monotonic function of t, which 
becomes obvious after one differentiates (p2[). 

Let us make estimations of the parameters during that 
intermediate stage. From (23) it is easy to find that w 
diminishes from 1 to a substantially smaller value during 
the time ~ 1/n. To have T/{Pe^w) of order unity 
by the beginning of the second stage, there should be 
w ^ l/inPe^). Looking at (|2^), we observe, that since 
H \s a. decreasing function of r, the left-hand side of ( |2^ ) 
is less than Hq w Sn? . On the other hand, the term 
2/w in the right-hand side can be estimated as nPe^. If 
n ^ Pe^, we can disregard the left-hand side. We as- 
sume in addition that the duration of the second stage is 
much less than r*. Then, we can substitute r by r* in 
the argument of and obtain the equation 



= 4 - 



12/1 



\Pe^wJ 



(25) 



We do not take square root, since w' changes sign dur- 
ing that part of the evolution so that both branches of 
( p5| ) are pertinent. Actually, (j2^) is valid for all r > r», 
since at w ^ r^/Pe^ it turns into ([2^). If our assump- 
tion about the duration of the second stage is correct, the 



transition region is described correctly too. On the other 
hand, if we add into ( p5| ) the term Hqw, which is small 
in the transition region, the equation thus obtained will 
correctly describe the evolution of w at all times r < . 

Therefore, we found that under the assumption of a 
short transition region, it is possible to reduce (|2|) to the 
equations w' = —^J%l^i{w) at r < and w' = v^'02(w) 
at r > T,, where 



^1=4- 
V'2 = 4- 



w 



Now we can continue with the second task, that is solv- 
ing (|l8|). Wc introduce also = w(t*) at the moment 
of reflection. Since at r = r* the derivative of w is zero, 
then we find from (G5h 



(26) 



Now (h8h can be rewritten in the following form 



2n 



w* V Pe^w^ 



dw 



dw 



'■VM^ ^^^^ 



(27) 



Estimating contributions into the integrals from the 
first and third time intervals, one can find that alone, 
they are too small to satisfy (p7|). On the other hand, it 
is easy to see, that if the derivative of ^/'i,2 at the point 
w* is not small enough, the second stage also does not 
contribute into the integrals in (|2^). The only way to 
have a solution is to make the derivatives small. Then 
"01,2 ~ V'l 2(^*)(^ ~ 'w*)^/2 in a wide interval, and both 
integrals in ( |27| ) logarithmically diverge, with cut off on 
non-zero value of the first derivative. Since we can make 
tjj'iwt:) arbitrary small by a small change of w*, the solu- 
tion exists. Equating the derivative of "0 to zero one finds, 
that is close to aT^,/Pe'^, where a is some number of 
the order unity, which depends on the form of (j) that is 
of the pumping x- The deviation of the first derivative 
from zero d is determined by (p7[): 



/ 2 1 

ij"iw,) S 



(28) 



One finds InS^^ oc n. Smallness of 5 justifyies our as- 
sumption on a short intermediate stage. Thus, we have 
shown that solution of the equation (18) exists, and 
y^l ^ l/{nPe^) which corresponds to ^^(O, 0) ^ iv'/^Pe. 
This gives the main contribution into (6'") oc ri'^"/^Pe" 
since the integral term in the action I is ~ n on our 
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instanton solution. Such moments correspond to the fol- 
lowing tail of the PDF oi lo = 0^ (note that it does not 
depend on the strain amplitude D) 



In [T{uj)] oc -{Ku^/Pf^. 



(29) 



Let us describe the above instanton solution in more 
physical terms. The instanton corresponds to some op- 
timal process giving the main contribution into (6'"). It 
produces large gradient which compensates for a small 
probability of such a process. Looking at (10,1^), we can 
distinguish the three stages of 9 evolution. During the 
first one, starting at t = — oo, the strain a ~ D stretches 
small-scale initial perturbation up to the width of order 
L and the force prepares some profile of 6 which has the 
amplitude of order y^P/D. The second stage starts when 
w is close to w*. Then a D,we can disregard both the 
advective and diffusive terms. Therefore, the width of 9 
does not change during that stage while the amplitude 
grows due to the force. Of all realizations of the force 
the constant one is preferred, since it gives the fastest 
growth. Then, 6 increases as 4>t. The weight of such a 
process is exp(— — t/2). The second term in the ex- 
ponent is the probability to have small a during the time 
t. Then we can find that 0^1 and t (x n (note that 
the second stage is long in terms of t yet it was short in 
terms of r). By the end of this stage 9^ oc n. And finally, 
during the last stage we can disregard the force. The 
profile having the amplitude n and width L by the be- 
ginning of the stage is compressed by the large negative 
a which can be estimated as tr ~ ~Dn. The duration of 
that stage (and the final width) is determined by diffu- 
sion: axdx9 ~ fid'^9 at the end. Then, the width of 9{x) 
is k/Dti while the amplitude is n^/P/D, therefore the 
final answer for (6*") is cx n^"/^(P/K)"/^, which corre- 
sponds to (p9[). To summarize, the optimal fiuctuation 



that gives the main contribution into (0") starts from 
an infinitesimal fluctuation which is initially stretched, 
then it has a long stage of suppressed advection when 
the amplitude of 9{x,t) grows yet it's spatial scale does 
not decrease, and then it is contracted fast. 

Let us stress that the instanton describes a very special 
configuration of the fields. For any other solution which 
does not satisfy correct self-consistency condition, there 
is no such a long intermediate stage of growth. If the 
parameters were not fitted to guarantee the existence of 
this long stage, either large a would bring a singularity 
in the solution or 9x would not be large enough at i = 0. 

Note that the scalar itself has an exponential PDF 
tail, the fact that the gradient PDF is less steep was 
correctly attributed in to the fluctuations of the dif- 
fusion scale. One may explain 2/3 stretched exponent 
in the following simple way: Due to diffusion, local gra- 
dient can be thought of as a product of a scalar fluc- 
tuation and inverse diffusion scale. The former has an 
exponential PDF tail [p| while the latter is proportional 



to the local stretching rate which is Gaussian. There- 
fore: (w") ~ (fp-y^l'^nP-r'^nP-l'^ which corresponds to ([2S|). 
Also, instanton formalism provides for an instructive in- 
sight into the relation between the scalar and it's gradient 
on the optimal fluctuation: comparing the second and 
the third terms in (^ one gets for the current dissipative 
scale oc ^JkI9D, substituting that into 9x — 9/rd we 
obtain Ox oc 9^/"^ so that exponential tail for 9 has to 
correspond to 2/3 stretched exponential for the gradient. 

To conclude, the gradient's PDF is given by (^ for 
Kup- jP < 1 and by (|^) for 1 <C ku? jP < Pe' which 
agree with the results found by a time-separation for- 
malism 1^. Speaking on generalizations, it is likely that 
one-dimensional instanton described here may be rele- 
vant for a multidimensional case, both compressible and 
incompressible, due to universality of a locally flat ramp- 
and-cliff structure discussed in [9|-|Tl|]. Note that the 
stretched exponential tail is what one expects for steady 
gradient's distribution (which is possible only when dif- 
fusion is present) [ p^ contrary to unsteady log-normal 
distribution which takes place without diffusion 0. 
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